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ABSTRACT 

We present high-resolution CO(l-O) observations of the lensed subminimeter galaxy (SMG) 
SMMJ14011+0252 at z = 2.6. Comparison to the previously-detected CO(3-2) line gives an in- 
tensity ratio of rs^i = 0.97 ± 0.16 in temperature units, larger than is typical for SMGs but within 
the range seen in the \ow-z ultraluminous infrared galaxy population. Combining our new data 
with previous mid- J CO observations, we perform a single-phase large velocity gradient (LVG) anal- 
ysis to constrain the physical conditions of the molecular gas. Acceptable models have significant 
degeneracies between parameters, even when we rule out all models that produce optically thin emis- 
sion, but we find that the bulk of the molecular gas has T^m = 20-60 K, riHs ^ 10^-10^ cm~^, and 
Nqo/ = lO^^-^^^^-^^ cm~^ km~^ s. For our best-fit models to self-consistently recover a typical CO- 
to-H2 abundance and a plausible degree of virialization, the local velocity gradient in the molecular 
gas must be substantially larger than its galaxy-wide average. This conclusion is consistent with a 
scenario in which SMMJ14011+0252 has a fairly face-on orientation and a molecular ISM composed 
of many unresolved clouds. Using previous Ha observations, we find that SMM J14011+0252 has 
a spatially resolved star formation rate vs. molecular gas surface density relation inconsistent with 
those of "normal" local star-forming galaxies, even if we adopt a local "disk-like" CO-to-H2 conversion 
factor as motivated by our LVG analysis. This discrepancy supports the inference of a star formation 
relation for high-z starbursts distinct from the local relation that is not solely due to differing choices 
of gas mass conversion factor. 

Subject headings: galaxies: high-redshift — galaxies: individual (SMM J14011+0252) — galaxies: ISM — 
galaxies: starburst — galaxies: star formation — ISM: molecules 



1. INTRODUCTION 



nerbauer et aI]|2QQ9| [Carilh et aL][2QTQ| [Aravena et aL] 



Since the discovery of a population of dusty, high- 
redshift submillimeter galaxies (SMGs), whose high star 
formation rates (SFRs) mean they contribute substan- 
tially t o the cosmic infrar ed background (see review arti 
cle by 



2010p . The SMG studies have shown that the CO (3-2) 
to CO (1-0) line ratio favors a value of r^^i ^ 0.6 (in 
units of brightness temperature) , indicating the presence 
of a multi-phase mole cular ISM that includes a subs tan- 



Blain et al 



20^1), significant efforts have been ^^l^ .^^^ "^^^7.^r^ ^Q^^J^^H^ 



made to char act erl z^TTfe sta r-forming m olecul ar gas [2010] [Iv ison et al.] pllj [ iJamelson et al. pOT Tp WEI^ 



(e.g., 
|et al. 



Prayer et a Lj|1998[ [T9 99| |Neri et aL 



2005; Tacconi et a l.||2006l |2008 



f. 



2003 
WHiTe 



Greve 



many 

of these studies have used solely mid- J CO lines (i.e., 
CO(3-2), CO(4-3), and CO(5-4)), recent expansions of 
the frequency coverage of radio telescopes into the Ka 
band have allowed observations of CO (1-0) — the best 
tracer of cold and low-density molecular gas — for SMGs 
whose redshifts fall near the z ~ 2 -3 peak of the distribu- 



tion for radio-identified samples ( [Chapman et al. 2005) . 



> 1 



poulations in other bands (e.g.. 


Papadopoulos et al.|2001 ; 


Carilli et al.||2002 Greve et al. 




K lamer et al. 
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similar line ratios have been observed in other types 
of high-r edshift star-forniing s ystems like Lyman break 
gal axies ([Riechers et al. 2010 ) and BzK-selected g alax- 
ies (Dannerbauer et al.||2009[ Aravena et al. 2010), and 
in low-redshift dusty galaxies (albei t with a much wider 
range of ra^i; e.g., Mauersberger et al. 1999[ Yao et al 



2003), this result is in direct contrast to those tor quasar 



host galaxies, whose near- unity val ue of r.^.i is consistent 



with s ingle-phase molecu lar gas (Riechers et al. 2006 
[2011a"; see also r^^i > 1 in jlvison et al.||2012p . 

While these CO (1-0) measurements have been essen- 
tial for determining the global multi-phase nature of the 
molecular gas in SMGs, the lack of high-resolution map- 
ping leaves some margin for error in interpreting line 
ratios. In order to be confident in assessing physical 
conditions in the gas, we need to know that the CO 
lines are being emitted from the same localized regions. 



2 



Sharon et al. 



Given the clumpy structure observed in S MGs (e.g., Tac- 
coni et al.||2008l iSwinbank et al.||2010), and the role 



major mergers are likely playing in the star formation 



bursts of (at least ) some of thes e syste ms (e.g.jConselice 
et al |[2QQ 3; Nar ayanan et al.|[2QlQ( |Engel et al. 2010 



Dave et al. 2010 ), it is not safe to assume a homoge- 
nous distribution of ga s conditions within every galaxy 



L] 120031 ). In addition, 
5MGs are gravitation- 



(e.g., |Mao et al.||2000l |Zhu et al 

since many of the well-studied 
ally lensed, it is possible that differential lensing (the 
variation in magnification factor across a spatially ex- 
tended so urce) could be affecting th e observed line ra- 
tios (e.g., Blain 1999 Serjeant] [2012]) . Indirect evidence 
for the spatial variation ol excitation in SMGs has been 
seen in the more extended distribution of CO (1-0) rela- 
tive to radio continuum emission (from the synchrotron 
radiation from supernova remnants, which trace m ore ac- 
tively star- forming and thus higher-excitation gas; |Ivison| 
et al [2011 ). The larger spatial extent of SMGs' cold gas 



reservoirs is also bolst ered by the larger line widths of 
CO(l -O) vs. CO(3-2) ( |Ivison et al.(2Qll| [Thomson et al. 
20T2|. 



The high-resolution CO observations necessary for 
evaluating spatial variation of CO excitation also en- 
able spatially resolve d com parisons between the Schmidt- 
Kennicutt relations (Schmidt 1959; [Kennicutt|i 998) seen 
at low and high redshilts. Determining galaxies' gas 
mass surface densities from CO surface brightness re- 
quires the assumption of a CO-to-H2 conversion factor. 
Observations of high, and sometimes unphysical, gas 
mass fractions in low- 2; IR-bright starbursts and some 



disk galaxy centers (e.g., iScoville et al.||1991 Garcia- 


Burillo et al.||1993l |Downes et al.| 1993| l^olomon et al. 


1997; Scoville et al. 1997( 


[Oka et al. 1998| Hinz & 


Rieke 2006; Meier et al.j20] 


LO) imply a range of aco, 



4.6M0 (Kkms-^ pc^) ; [Solomon fc Barrett' l99l| ) for 
normal star-forming galaxies and a '' star burst" value 

(acq = O.8M0(Kkms" 



|d^ 



Downes & Solomon 



1998) for mergers and other high surface density envi- 



ronments. Some effort has be en made to find a suitable 
continuous form of ac q (e.g.. Glover & Mac Low"2011^ 
Narayanan et al.|2012 ) that also matches obs erved trends 



in inetallicity (e.g., [W ilson 1995; A rimoto e t al. 1996; Is 



ra^[T997l |Leroy et^.|,2006, .Bolatto et al.||2008, .Genzel 
et al.|| 2012| ). Which value (or values) of aco is appro- 



priate lor high-redshift galaxies is debated (e.g., Tacconi 



appro- 



|et al. 2008), and the apparent bimodality in the high- 2; 
Schmidt-Kennicutt relation (where SMGs fall along the 
\ow-z starburst track, and more "normal" high-z galax- 
ies extend the low- 2; disk galax y track to higher masses 
and star forma tion rates; e.g., Daddi et al.|[2010| [Gen- 
zel et al.||2010l cf. jlvison et al 



201ip may be caused 



in part by the choice ol aco- Comparisons between the 
Schmidt-Kennicutt relations at high and low redshifts 
may also be affec ted by the d ifferences in the CO lines 
observed (Narayanan et al.j2011 ), as the molecular gas at 
low redshilt is measured via the CO (1-0) line, while the 
molecular gas at high redshift has typically been mea- 
sured via mid-J CO lines. While galaxy-wide measure- 
ments of the star formation rate and molecular gas mass 
are typical for high-z systems, evaluation of the true sur- 
face density version of the Schmidt-Kennicutt relation 



(e.g., Kennicutt et al.j 


2007 


Bigiel et al.[[2008 Wei et al. 


2010[) requires spatial 


ly resolved maps ol CO lines and 



In order to explore the spatial variation of CO line ra- 
tios and the Schmidt-Kennicutt relation for high-z star- 
bursts, we have obtained high-resolution Karl G. Jansky 
Very Large Array (VLA) mapping of the CO (1-0) line 
in SMM J1401 1+0252 (J14011 hereafter), which we com- 
pare to previous mid-J an d rest-frame mid-IR observa- 
tions. J14011 (zopt = 2.55; Barger et al.[[199 9') was dis; 



covered as part of the SCUBA Lens JSurvey (Small et al.[ 
|1998[ [2002); it consists of two optical/near-IR compo- 
nents m the source plane (Jl and J2) that are weakly 
gravitationally lensed by th e cluster A1835 (z = 0.25) 
and one of its members (Jlc; [Motohara et al.|2005 ) , with 
a median magnification factor ol /i ^ 3.5 ± 0.5 ( [Smail[ 
et al.[ [2005). Since its discovery, J14011 has been fol- 



lowed up extensively w ith optic al and near-IR observa- 
tions dBarger et al.[[l999 ; Iviso n et al.|[2000[ [20011 [Tecza 
et al.|2UU4||I^rayer et al.i2004uSwiirbank et al.|2U04[ |Mo- 
tohara et al.||2005p Nesvadba et al. [2007^ , and it was 
only the second SMG to be successlully detected in any 
CO transition (CO (3-2); Frayer et al. 1999). Further 



CO(3-2) observations were made by Ivi son et al. (|200l| 
as were CO (3-2) and CO (7-6) observations by [Downef 



^and CO(7-6) observations by [Uownes 

fc Solomon] ([2003| . Strong upper limits on its X-r ay ILiix 
( [Ivison et al. [[ 2000 [ [ Fabian et al [[2000| [Rigby et al.](2008 ) 
and ratios ol the 7.7 fim polycyclic ar omat ic hydro car b on 
(PAH) feature to mid-infrared flux (Rigby et al. 20Q8| 
indicate that J 140 11 is a pure starburst, lacking a sigmi- 
icant contribution to its dust luminosity from an active 
galactic nucleus ( AGN) . It is si milar in its far-infr ared lu- 
minosity (I/piR ~ 4 X 10^^ L0;|Smail et al.[[2005[) and its 
clumpy morphology within and between optical/near -IR 
components (Ivison et al.n2000j , Tecza et al. 2004, Swin-[ 
l.|2Q04p to low-redshitt ultra/ luminous inlrared 
galaxies ( U /L iRGs) su ch as Arp 220 d^rayer et al.[|1999l 
Ivison et a"L][2Q 00; Do wnes &: !Solomon[|2003[ IMotohara 
' IRigby et al.[[2008p. 

hrcD]i-o) liAj-mai- 



et al.l 



J2005] ]Nesvadba et a 

Initial single-dish observa tions ol t 

cated rs^i = 0.76 ± 0.12 ( Harris et al.||2010 [), similar to 
the values seen for other SMGs. 

We assume the WMAP7+BAO+i^o mean ACDM cos- 
mology throughout this paper, with = .725 and 
i^o = 70.2kms-iMpc-^ ^Komatsu et al.||201l[ ). 



2007 



2. OBSERVATIONS 

We observed J14011 (a(J2000) = 
^(J2000) = +02°52^24 .1^^ CO(3-2) position from 



Downes fc Solomon |2003[ ) at the VLA on 2010 October 
1 and 3 in the DnC configuration (maximum baseline 
1.868 km) and 2011 January 22, 25, and 28-31 in the 
CnB configuration (maximum baseline 6.903 km). Dur- 
ing the DnC observations, 23 antennas were functioning 
and equipped with Ka-band receivers, and the weather 
conditions were mixed (< 70% sky coverage by cumuli- 
form clouds). During the CnB observations, 26 Ka-band- 
equipped antennas were available, and the weather con- 
ditions were excellent. We observed with the Wideband 
Interferometric Digital Architecture (WIDAR) correlator 
in the "Open Shared Risk Observing (OSRO) Full Polar- 
ization" mode using the lowest available spectral resolu- 
tion (64 channels x 2 MHz resolution). The intermediate 
frequency (IF) channel pairs A/C and B/D were tuned 
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3 



to overlap by ten channels, creating a full bandwidth of 
236 MHz centered at 32.3703 GHz (chosen to keep the 
line away from IF pair edges). At the beginning of each 
track we observed 3C286 as a flux calibrator, adopting 
Sy = 1.9455 Jy as the Common Astronomy Software 
Applications (CASA) package's default "Perley-Butler 
2010" flux densit}F] Phase and amplitude fluctuations 
were tracked by alternating between the source and a 
nearby quasar, J 1354-0206, with a cycle time of 6 min- 
utes; this quasar was also used for passband calibration. 
A total of 18.7 hours was spent on source. 

Calibration and mapping were carried out in CASAH. 
After the tracks had been corrected for Doppler shifts 
relative to the local kinematic (radio) standard of rest 
using the CASA routine cvel, the original channels 
were smoothed by a factor of two, resulting in 4 MHz 
(37kms~^ rest frame) spectral resolution. The natu- 
rally weighted channel maps have a synthesized beam 
of 0.^^68 X 0.^^52 at position angle 119.02 degrees and an 
average RMS noise of 53.2 /iJy beam~^. Analysis of the 
resulting data cube was performed using a custom set of 
IDL scripts. 

3. RESULTS 

Our spectrum (Fig.jl]) successfully detects the CO (1-0) 
line at the 8<j level. A Gaussian fit to our line pro- 
file gives a peak flux density of 2.08 ± 0.22 mJy and 
a FWHM of 151 ± 19kms-^ The FWHM is consis- 
tent with the CO(7-6) line width of 170 ± 30kms-\ 
and marginally narro wer than the CO (3-2) line w idth of 



190 ± 11 kms~^ from [Downes fc Solomon] ( |2003[ ). The 
integrated CO (1-0) line tlux (using the seven channels 
around the line peak) is 0.32 ± 0.03(±0.03) Jy kms~^, 
where the latter 10% uncertainty is associated with the 
flux calibration. While our measured integrated line flux 
appears to be lowe r than the ^ 0.40 ± 0.05 Jykms"^ of 
Harris et al. (12010), the two values are consistent within 
their la errors. The slight discrepancy is likely because 
the Zpectrometer-measured value is from a Gaussian fit 
in which the line width is only an upper limit (the nar- 
row line was marginally spectrally resolved by the Zpec- 
trometer). Convolution of the VLA spectrum with the 
20 MHz FWHM sinc(x) response function of the Zpec- 
trometer reproduces th e Gaussian fit integrated line flux 



given in Harris et al.| ( [2010^ , indicating that our mea- 



surements are m good agreement when we account for 
spectral resolution differences. Using the Milky Way con- 
version factor ag o = 4.6 M(t^ (Kkms~^ pc^) ^(Solomon 
fc Barrett 1991) suggested by our LVG analysis (Sec 



tion^TT) and correcting for lensing magnification and he- 
lium, the molecular gas mass is (1.9±0.3) x 10^^ Mq. The 
redshift as determined by our CO (1-0) measurement is 
z = 2.5653^0.000 1, which is consistent with the previous 
optical (e.g., Bar ger et al. 1999; Ivison et al. 2000) and 
CO (e.g., jPownes !Solomon,,2003, ^Harris et al.„2010^ 
redshift s. 

The integrated intensity map (Fig.[2| shows that the 
CO (1-0) emission is spatially extended in a north- 
east/southwest direction, parallel to the lensing shear 
( [Smail et al.^^2005) , and contains two peaks embedded 
within more extended emission, similar to those seen 
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Fig. 1. — STkms"-*^ resolution VLA spectrum plotted rela- 
tive to the z = 2.5652 CO(l-O) systemic redshift from the 
GBT/Zpectrometer observations ( ^Harris et al.|2010| ). Vertical lines 
indicate statistical uncertainties. 
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Fig. 2. — Integrated CO(l-O) intensity map. Contours are mul- 
tiples of ±2cr (cr = 0.16 mJy beam" ■'^), where negative contours are 
dotted; synthesized beam is shown at lower left. 



in Ra (Tecza et al. 2004). The center of the CO(l-O) 



line emission is consistent with the peak positions of 



the CO (3- 2) and C O/7-6) hnes 



Downes & Solomon] ([2003|) (Fig.[^ 
cal 



reviously detected by 
] with cTi 0.''05 
t in the image plane 



based on an elliptical Gaussian 
and the intrinsic uncertainty of the phase calibrator po- 
sition. The Gaussian fit gives a major axis FWHM of 
2/'6 ± O.^'l at position angle 18.7 ± 1.3 degrees and a 
minor axis FWHM of 0'^8 ± 0.'''1 when accounting for 
the beam size. As J14011 is lensed, these dimensions do 
not easily translate into rest-frame physical dimensions. 
While the line FWHM is quite narrow, the CO(1-0) emis- 
sion does show a velocity gradient. In Fig.[4j we plot 
overlaid channel maps of the red and blue halves of the 
central 200kms~^ of CO(l-O) line as well as the same 



in 


Tecza et al. 


(2004). 


lines are slightJ 


y ottset 



^ http://casa.nrao.edu 



peaks in the integrated line maps apparently dominated 
by emission at different velocities. In Fig.|4j we also show 
a CO(l-O) position-velocity diagram for J14011 for a po- 
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Fig. 3. — Overlaid contours of all published CO maps of J14011 
and the Aguirre et al.l ( |2012| > HST F160W image (with the 
foreground mterloper, Jlc, removed), centered at q;(J2000) = 
;L4hoimo4.sg3 (5(J2000) = +02°52'24/'l. The CO(l-O), 

CO (3-2), and CO (7-6) maps are in yellow, blue, and red respec- 
tively; the synthesized beam for each map is shown by the filled 
ellipse of the corresponding color in the lower left corner. The 
CO (1-0) line contours are as described in Fig.jg] The COf3- 
2) and CO (7-6) observations are taken from |Downes Solomon] 
(|2003|. The CO (3-2) contours are multiples ot ±2(7, where a = 
0.41 mJy beam" The CO(7-6) contours are in steps of ±lcr, 
starting at zb2cr, where a = 1.3 mjy beam" -"^ . 

sition angle of 75 ± 5 degrees that maximizes the velocity 
gradient. 

We do not detect any significant continuum emission 
from J14011 in the (rest-frame) 2.6 mm continuum map 
produced from the line- free channels. We estimate a 3cr 
upper limit for the continuum emission of 26.7 /iJy for a 
point-like source, or 157 /iJy when scaled to the dimen- 
sions of the CO (1-0) emission. 

In order to compare the CO (1-0) emission to existing 
optical and infrared data, we have taken care in aligning 
images. Based on the source sizes and structures in the 
two sets of maps in Fig.|4| and the fact that Ha emis- 
sion from Hll regions is expected to follow molecular gas 
at this linear resolution, we can use the positions of the 
northern and southern peaks to better align the datasets. 
The average offset between the peaks was applied as a 
shift to align the SINFONI Ha data cube with the radio 
observations. The measured offset is 0.^^82, which cor- 
responds to roughly ^ 2x the SINFONI point spread 
function (the astrometric uncertainty of the CO (1-0) 
map is ~ 0.^^09, including both statistical uncertainty 
and uncertainty due to the phase calibrator position). 
Based on this alignment, we use the continuum emission 
from the line- free channels of the SINFONI data-cube (at 
- 2.2 /im) to ex tend the alignment to the HST/mCMOS 
F160W data of [Aguirre etal.| ( |2012[ ). As we expect the 
continuum emission at both wavelengths to come from 
the same physical regions in the source, we use the con- 
tinuum from the foreground interloper (Jlc) and the UV- 
bright companion galaxy (J2) to determine the shift to 
apply to the iJST image to align it with the Ha and thus 
CO(l-O) maps (Fig.[3|. We measure an offset of 0.^^49 
between the F160W and 2.2 /im maps, making the offset 



between the CO(l-O) and iJST maps < 0.''5. 

4. ANALYSIS 

4.1. LVG modeling 

J 140 11 has integrated brightness temperature ratios 
r3,i = 0.97 ± 0.16 and r7,i = 0.20 ± 0.04. Comparisons 
between the higher-J and CO(l-O) emission (convolved 
to the resolutions of the higher-J maps) do not show ev- 
idence for spatial variation of CO excitation. This value 
of rs^i is larger than the typical rs^i ~ 0.6 for SMGs 
and more similar to the r^^i ~ 1 seen in high-z quasar 
host galaxies. The strong upper limits on X-ray emis- 
sion and observations of mid-IR PAH features indicate 
that J 140 11 likely does not contain an AGN of sufficient 
strength to influence the global gas conditions, making 
its near-unity r^^i unusual for its class. However, among 
low-z U/LIRGs (which also have an average r^^i ^ 0.6), 
there are a small number of starbursts that globally or 
locally have r.s.i 1 but lack AGN, such as N GC 3690 
dlono et al.||20Q9| ), UGC 873 9 (|Yao et al.| |2QQ3|, the in- 
ner region olNUC 253 (e.g. 'Bayet et al."'2004V, and the 
northeast lobe of M82 (e.g.. Ward et al. 2003). While the 
high r^ A of J14011 is ra re for a pure star burst, it is not 
unique. Yao et al. ( 2003) propose that high CO excita- 



tion in U /LIKGs can be caused by a very high conce ntra- 
tion of star f ormation, consistent with the finding of'I ono| 



et al. (|2009|) that the peak central rs^i is greater than or 
equal to the global rs^i in spatially resolved maps. If 
the duration of highly concentrated star formation (such 
that the global and "peak" CO excitation are identical) 
is short-lived, then it is sensible that only a small number 
of pure-starburst SMGs (or local U/LIRGs) will exhibit 
high ra^i. Gas distributions that are compact relative to 
heating sources may be the trait that such systems share 
with quasar host galaxies. 

In order to investigate the physical conditions of the 
gas that produces these ratios, we have undertaken an ex- 
citation analysis using a Large Velocity Gradient (LVG) 
model. The details of the LVG modeling and analysis are 
described in Sharon et al. (in prep.). In short, the model 
follows the method of [Ward p002) , assumes a spheri- 
cal cloud geometry, ana is moditied to include the effects 
of the cosmic microwave background (CMB). With ob- 
served CO transitions equal in number to the model pa- 
rameters (kinetic temperature, H2 density, and CO col- 
umn density per unity velocity gradient) and substan- 
tial uncertainties in the measured line ratios, we expect 
significant degeneracies between the allowed model pa- 
rameters. We therefore use both a minimum ^ind a 
Bayesian a nalysis (the latter better capturing our un- 
certainties; IWardetaL 2003) to determine which model 
parameters best reproduce our measured line ratios. For 
the Bayesian analysis, we have applied logarithmic priors 
to the LVG model parameters to compensate for the large 
logarithmically-sampled parameter space, as opposed to 
the usual agnostic prior equal to unity. 
The 



parameters of the minimum x model are 
136 K, = 10^-2 cm-3, and Nco/Av = 

10^^"^^ cm~^ km~^ s. However, there are significant de- 
generacies in all three parameters. For example, the 
Bayesian maximum probability model favors the oppo- 
site extreme in temperature-density parameter space: 
T = 24 K, = 10^-^ cm-^ and Nco/Av = 
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Fig. 4. — Left/Center: Overlaid contours of two 10 MHz (~ 110 km s"-*^) channels centered at 56.52 km s"-*^ (red lines) and —54.75 km s"-*^ 
(blue lines) for the CO(l-O) line (left) and the Ha line (center). Contours are multiples of ±2cr (cri-o = O.lTmJy beam"-*^; (jHa = 
1.18 X 10"-*^^ W/im"-*^ m~^), where negative contours are dotted; synthesized beam/point spread function is shown at lower left. Right: 
Position-velocity diagram of the integrated line where the major axis (along the velocity gradient) is at a position angle of 75 degrees that 
maximizes the velocity gradient. Contours are multiples of 1.75 mJy beam" -"^ . 



-|16.75 , 



produced for 20 K 



columns per velocity gradient (10 



Similarly well-fitting models are 
5, T S, 150 K at comparable CO 



17±0.25 



cm ^ km -'^s), 



where lower temperature models favor higher H2 densi- 
ties 10^ cm~^) and higher temperature models favor 



lower H2 densities (~ 10 cm 



The optical depths 



produced by these two models are very different from 
one another. The minimum model produces optically 
thick emission for the CO (3-2) and CO (7-6) transitions 
(t3 2 = 13.4, Tr 6 = 11-4), and optically thin emission 
for the CO(l-O) transition (ri = 0.7). The cooler and 
more dense maximum probability model has the optical 
depths of the CO (1-0) and CO (7-6) lines reversed with 
Ti = 2.0, Ts 2 = 11-8, and tj q = 0.7. We also consider 
a restricted parameter space for the LVG models in which 
all of the measured lines must be produced by optically 
thick emission. The minimum model with Tj^pper<7 > 
1 has a significantly lower temperature, T = 60K, = 
10^ '^ cm-^, and Nqo/Av = 10^^-^^ cm-^ km"^ s, while 
the maximum probability model only slightly increases 
in temperature (relative to the previous maximum prob- 



28 K, 

-1 . 



10 



5.0 , 



and 



ability model), with ^ 

Nco/Av = 10-^^-^^ cm~^ km~^ s. In these models, the 
optical depths of the CO (1-0) and CO (3-2) emission are 
comparable (ri ^1.5, rs 2 ^ 10-15), but the CO (7-6) 
line is only just optically thick in the colder model, while 
in the warmer model the optical depth is a factor of 
~ 10 larger. The CO spectral line energy distributions 
(SLEDs) for these four models are shown in Fig. [5] 

Our cooler and optically thick models favor ki- 
netic tempera tures comparabl e to the dust temperature, 
Tdnst = 42K (iWu et al."2009), and are in li ne with the 



results of other LVG analyses of S MGs (e.g.,|Weifi et aL 
20071 ICariUi et al.|[2Q10t Riechers et al. 2010; Daniel- 
son et al. [20111 IRie chers et al. 2011b) . Observations 
of ditterent molecular or atomic species can be useful 
for b reaking temper a ture d egeneracies in CO LVG mod- 



els. [Walter et al. (2011) make an independent esti- 



mate of the kinetic temperature in J14011 using their 
C i(^P2 ^Pi) line measurement in conjunction with 
the p reviously-measured C i(^Pi ^Po) line ([Weifi et al. 
^QQ5t , deriving Tci = 32.4 ± 5.2 K. At fixed TkhT^^TcTT 
the minimum model (which is indistinguishable from 
the maximum probability model due to similar probabil- 
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Fig. 5. — Measured line ratios (points; in flux units) and CO 
SLEDs for the maximum probability model (thick solid line), max- 
imum probability model requiring r > 1 (thin solid line), minimum 
model (thick dashed line), minimum model with r > 1 (thin 
dashed line), and the maximum probability /minimum model 
at flxed Tci (thick dotted line). The LVG model parameters that 
produce these SLEDs are as labeled. 

ity and distributions in the nY^^-Ncol lS.v plane) has 

1 o 



= 10^-^ cm ^ and Nqo/^v 



(Fig.[5| and produces optically thick emission for all three 
transitions (ri o = 2.1, ra 2 = 16.0, r7_6 = 2.7). While 



a prior LVG analysis performed by Bayet et al. (2009) 
preferred higher temperatures (T = 205 Kj and lower 
densities {nn^ = lO^-^cm"^) than any of our well-fitting 
models, due to the different goals of their work (pre- 
dicting line strengths for a statistical sample rather than 
accurate determination of individual galaxy properties), 
given their use of only the CO (3-2) and CO (7-6) line 
measurements and our consistency with other estimates 
of the gas temperature, we view our results as more ro- 
bust. 

Besides the degeneracies between the best-fit models 
noted above, it is worth considering what degeneracies 
exist within each model (e.g., degeneracies between 
and Nco/ Av for a specific choice of temperature, or be- 
tween Tkin and for a specific value of Nqo/Av). In 
Fig. [6] we have plotted curves of constant probability (or 
X^) on the temperature-density plane for the best-fit val- 
ues of Nco/Av. For a specific value of Nqo/Av, there 
is generally an "L" -shaped region with nearly constant 
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probability /x^ that will reasonably reproduce a given 
CO line ratio. With multiple line ratios, the locus where 
the two "L" -shaped regions overlap produces the best- 
fit models. For multiple line ratios, these "L" -shaped 
regions overlap more for emission lines with similar op- 
tical depth ratios. The r > 1 regions are marked with 
dashed lines in Fig. [6] and are also "L" -shaped. As optical 
depth is proportional to Nqo/Av, larger regions satisfy 
the r > 1 criterion for larger values of Nqo/Av. 

For a specific temperature, there are also degenera- 
cies between the H2 density and CO column density per 
velocity gradient. In Fig. [7] we have plotted curves of 
constant the nu^-Nco/^v plane for the optically 

thick model temperatures that give highest probability 
and minimum (i^i this plane, the probability and 
contours have the same shape, but are scaled differently 
since the probability is ex ). While slightly shifted 
between the panels for the two temperatures, there is a 
clear linear valley with low x^- This region corresponds 
to Nco/^v oc ^H^' exists for all Tkin > 30K 

models (for T < 30 K, the valley of low values asymp- 
totes to a single value of Nqo/^v for higher H2 densi- 
ties). This linear feature corresponds to a constant opti- 
cal depth ratio for the emission lines that best reproduce 
the measured line ratios. 

With these substantial degeneracies it is worth con- 
sidering priors in addition to optical depth that might 
limit the allowed parameter ranges for the LVG mod- 
els. Since the most dense molecular gas in the Galaxy is 
found in star-forming virialized clouds, we can consider 
placing limits on the density by considering the degree 
of virialization as parameterized by 



1.54 



[CO/H2 



CO 



10^ cm~ 



-1/2 



(1) 



where Aco is the CO abundance per unit velocity gradi- 
ent and a ~ 1-2 .5 is a constant that depends on the c loud 



geometry (e.g., |Goldsmith||2QQl| |Greve et al.||2QQ9D . In 

this parameterization, [CO/H2] A^q is the velocity gra- 
dient in the molecular gas, and K^ir <C 1 is dynamically 
unobtainable. For J14011, the velocity gradient could 
be a value determined by global parameters (given by 
the - 150 km s-^ CO(l-O) FWHM and app roximate de- 
lensed optical source size of ~ lOkpc from | Small et al.| 
([2005)) if the line luminosity is trac ing the total gravita - 
tional potential of the galaxy (e.g., Downes et al.||1993). 



or a different value if the molecular gas is in individual 
star- forming clouds. For the global velocity gradient, re- 
quiring that Kvir > 0.1 means that nn^ > 10^-^ cm~^ is 
not allowed, excluding all of the well-fitting models dis- 
cussed above. In general, low x^ models with low den- 
sities that meet the i^vir > 0.1 criterion under-predict 
the CO(3-2)/CO(l-0) line ratio by > lcr3,i. Velocity 
gradients larger than the global value would allow the 
higher densities of our best-fit LVG models to satisfy 
the virialization requirement. In order to allow densi- 
ties as high as 10^ cm"^ (10^ cm"^) with i^vir > 0.1, the 
velocity gradient would need to be > 0.27kms~^ pc~^ 
(> 0.86kms~^ pc~^.) If the individual clouds are virial- 
ized (i^vir = 1)7 these H2 densities require velocity gra- 
dients of ~ 3kms~^ pc~^ (~ 9kms~^ pc~^). 
Models with low uyl^ and high Nqo/Av can be elimi- 



nated by requiring the CO column length to be less than 
the maximum source size. However, this prior requires 
the assumption of an upper limit for [ CO/H2]. Conser- 
vative limits (i.e., [CO/H2] < 5 x lO-^^^Ward et al.|2003|) 
exclude clouds with densities comparable to or less than 
that seen in the Milky Way (~ 10^ cm~^) along the low- 
X^ degeneracy valley in the nu^-Nco/Av plane. While 
this prior is not useful for distinguishing among the best- 
fit models for J14011, it does support the higher H2 den- 
sities preferred by the LVG analysis. 

While it is unlikely that a simple LVG model will 
accurately capture the full range of temperatures and 
densities that exist in real molecular clouds, consider- 
ing the degeneracies discussed above, we conclude that 
the bulk of the molecular gas likely has, on average, 
Tkin = 20-60 K, ^ lO'^-lO^ cm-^ and Nco/Av = 
-|^q17.oo±o.25 i^^-i g rpj^^ maximum likelihood for 
the thermal pressure peaks at ~ lO^Kcm"^, but has 
a significant tail up to lO^Kcm"^. Measurements of 
other mid- J CO lines, like CO(5-4), would help to 
better constrain these models. Using our best esti- 
mates for the LVG model parameters, we constrain 
[C0/H2]/(Av/Ar) = lO-^ MO-^-^pckm-^s, which is 
close to the standard assumption of 10~^ pckm~^ s used 
in other LVG analyses. While the measured line ratios 
in J14011 appear consistent with a single phase molecu- 
lar ISM, analysis of other SMGs (e g., [Harris et al. 2 QTq| 
Ivison et al1[2QTT| |Danielson et al.||2Ull| lEothweh et^ 
201 2D and observations of local galaxies (e.g., lYao et al. 



2003) indicate that gas in star-forming galaxies is multl 



phase; i.e., it is unlikely that the gas responsible for the 
bulk of the CO (1-0) emission is the same gas that is 
emitting the CO (7-6) line, which is effectively what we 
are requiring by considering only single-phase models. 
With only three measured lines, the results of a multi- 
phase LVG analysis would be far too degenerate to draw 
meaningful conclusions about the state of molecular gas. 
Higher resolution observations would be helpful in estab- 
lishing whether the CO lines are or are not being emit- 
ted from the same physical volumes of gas. Based on the 
agreement of the global line profiles for the existing CO 
measurements it is unlikely that excitation differences 
exist on large scales across J14011, but local variations 
with broadly similar kinematics could exist on scales be- 
low our current resolution limit. 

Based on the CO abundance per unit velocity gradient 
of the best-fit LVG models, we can determine the CO-to- 
H2 abundance for an assumed velocity gradient. Using 
the global velocity gradient from our measured line width 
and the Small et al. (2005 ) source size, we estimate the 
CO-to-H2 abundance to be [CO/H2] = IQ-^-^-K)-^-^ , 
which is very low. Alternatively, if we assume more 
standard values of [CO/H2] = lO-^-lO""^, the best-fit 
model densities return higher velocity gradients of ~ 1- 
30kms~^ pc~^, a range that includes the necessary gra- 
dients for our best-fit models to meet the virialization 
criteria. 

Since LVG modeling produces reasonable physical con- 
ditions for J14011, and it is only the introduction of the 
global velocity gradient that produces very low CO-to-H2 
abundance and unphysical virialization parameters, we 
conclude that J14011 likely has large local velocity gra- 
dients. Higher velocity gradients are possible if J 140 11 



CO(l-O) in J14011 



7 





Fig. 6. — Contours of constant probability (left) and (right) in the temperature-density plane for the best-fit values of Nco/Av. 
Models with higher probability than the highest contour (or lower than the lowest contour) have been logarithmically shaded in 
greyscale. Models with parameters outside of the green-shaded region produce optically thin emission in at least one line. The best-fit 
models with and without the optical depth cut have been marked with yellow crosses and circles respectively. 
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Fig. 7. — Contours of constant i^i the nH2~^Co/^'^ plane for highest probability optically thick model (28 K; left) and minimum 
optically thick model (60 K; right). Models with < 0.1 have been logarithmically shaded in greyscale. Models with parameters that fall 
below the green-shaded region produce optically thin emission. The highest probability /minimum models in each plane are marked by 
the yellow crosses. 

has a face-on orientation and contains numerous small 
clouds below our resolution limit that have turbulent ve- 
locities and/or bulk motions summing to the measured 
FWHMi ^ 150kms~^. If the individual clouds are 
virialized, their velocity gradients are large enough to 
produce more reasonable values of the CO-to-H2 abun- 
dance for our best-fit densities. Although such velocity 
gradients would cause J14011's molecular clouds to fall 
above the Galactic line-width vs. size relation ( [Larson 
11981'), such behavior has been seen in another 2; ^ 2 
Wg, SMM J2135-0102, which has a mean internal ve- 



l ocity gradient of ^ 0.5 k m s ^ pc ^ for individual cl umps 



(Swinbank et al. 2011). Swinbank et al. (2011) also 



order to convert the Ra surface luminosity to a star for- 
mation rate surfac e density, we use the Lhq!-SFR conver- 
sion factor giv en in|Kennicutt| (p^998 ) scaled to match the 
'Calzetti et al."(^007) initial mass function (IMF) used in 
Bigiel et al. (2008, 2010). We have also applied an ad- 
ditional correction to reflect the fact that the observed 
Ra luminosity underestimates the galaxy's total star for- 
mation rate: base d on the [K ennicutt (1998) calibration 
rescaled to the C alzetti et al.| ( |2007| ) IMF, J14011's to- 
tal FIR-determined SFR is 4OOM0yr~^, exceeding its 
Ha-determined SFR by a factor of fifteen. Given that 
J 140 11 does not contain a strong AGN (Rigby et al 



find similarly high densities, and thus infer large velocity 
gradients under the assumption of virialization, for the 
molecular gas in SMM J2135-0102. 

4.2. Resolved star formation law 

Using the aligned Ha and CO (1-0) maps, we plot a 
spati a lly resolved Schmidt-Kennicutt relation (Schmidt 



2008|) , that there is good agreement between the CO and 
the Ha morphologies (unlikely if the FIR vs. Ha discrep- 
ancy were determined by large-scale patchiness in obscu- 
ration), a nd that the c orre ction for the Balmer decre- 
ment as in Tecza et al. ([200 4 ) would already account for 
lie facte 



1959] |Kennicutt| |1998) for J14011 in Fig.js] We have 
smoothed the integrated line maps to the sai 



maps to tJie same spatial 
resolution (0.'''70 x 0.^^65), resampled to the pixel size 
from the SINFONI data (0.''25), and excluded pixels that 
do not have at least 2a significance in both maps. In 



a factor of 1.5 (of the factor of 15 difference), we have 
simply scaled the Ha surface brightness of each pixel by 
the ^/o6a/ ratio SFRfir/SFRhq; before plotting in Fig.js] 
Converting CO (1-0) luminosity to molecular gas mass 
requires the assumption of a conversion factor aco- 
Based on the results of our LVG analysis above, which 
produces more self-consistent molecular gas properties 
if the molecular ISM in J 140 11 is composed of mul- 
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Fig. 8. — Star formation rate surface density as measured by Ho; surface brightness (rescaled by the glob al SFRftb /SFRnrv ) vs. molecular 
gas mass surface density for J 140 11 (crosses) and total gas mass surface density for local disk galaxies from Big iel et al.| ( |2010| ) (dots; clipped 
at 2 X 10~^ Mq yr"-*^ kpc~^). The average value for J14011 is given by the star, and a power law fit to the resolved surface density points 
is given by the solid line. The arrow pointing to the left shows where J14011's locus of points (centered at the average value) would move 
if the molecula r gas m ass were calculated using the "starburst" value of aco instead of the "disk" value. Dotted lines are the same as 
in [Bigiel et al.| (12008*), where the diagonal lines represent the constant star formation efficiencies required to consume 1% (yellow), 10% 
(orange), and 100% (red) of the gas in 10^ years, the left vertical line shows the Hi surface density saturation, and the right vertical line 
marks the proposed transition between star formation laws at 200 pc~^. 



tiple unresolved clouds rather than a single large gas 
reservoir, we use a fiducial Milky Way disk value of 

acn 4.6M0 (Kkms"Vc^)"^ ( jSolomon fc Barrett 



1991). This choice is consiste nt with the propo sed con- 



tmuous form of aco of Narayanan et al. ( 2Q12[) for the 
measured metallicity of component Jl froml'Tecza et al.' 
(^2004) of 12 + log(0/H) = 8.961^;?^. In Fig.g the 
gas mass surface densities for both J 140 11 and the low- 
redshift galaxies have been scaled to account for helium. 
Since gravitational lensing conserves surface brightness, 
magnification corrections are to first order unnecessary 
for evaluating the resolved Schmidt-Kennicutt relation 
for J14011. 

The resolved star formation law for J14011, as probed 
by the FIR-corrected Ha map, is not consistent w ith that 
of the local disk galaxies from "Bigiel et al. (2008) despite 
our use of a Milky W ay gas mass conversion factor. While 
Bigiel et"ar| ( |2008 ) use a sub-Milky Way value of aco = 

3.2M0 (Kkms~^pc^) ^, any value of aco < 4.6 moves 
J 140 11 further above the local relation, and aco =0.8 
m oves it far above. Ou r result supports the conclusions 
of Genzel et al. (2010), indicating that SMGs' position 
above the local star formation relation is not simply due 



to the choice of CO-to-H2 conversion factor. Correcting 
for different IMFs while respecting the differing choices in 
aco, the star formation law for J 140 11 is consistent with 
both the SMGs in Genzel et al.| (|2010|) (who use aco = 

( |1998D (who 
the high-z 



1.0) and the local star bursts froni ]Reniiicutt 

uses aco = 0-8), and is not consist ent wit 

"normal" star-forming galaxies from| Daddi et aL ]( 2010[ ) 
(who use aco = 3.6). Re scaling J14011 to the lower 
aco used in Genzel et al. (2010) or Kennicutt^ (|1998[) 
would move J14011 far above the star formation relation 
for SMGs and low- 2: starbursts, further supporting our 
choice of a Galactic CO-to-H2 conversion factor for this 
system. 

To allow comparisons to other work, we fit the star 
formation and molecular gas surface densities to a power 
law and formally retrieve a near-unity index of 0.98 ±0.1 4 
(albeit with = 2.48), consisten t with indices measured 
in local spirals ( Bigiel et al.|2008' ). Our application of the 
scaling factor to account for dust-obscured star formation 
does not affect our measured index, as appropriate if Ha 
traces the location of J14011's star-forming regions, even 
if it does not capture the total SFR. We note that in 
our pixel-to-pixel comparison, the data points are not all 
independent due to our oversampling of the beam/point 
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spread function; the points in Fig.jS] are drawn from a 
total area of ^ 10 resolution elements. 



5. CONCLUSIONS 

We present high-resolution CO (1-0) mapping of the 
submillimeter galaxy SMMJ14011+0252. Based on a 
spectrally-resolved detection, we find a revised near- 
unity value for the CO(3-2)/ CO(l-O) line ratio (rs^i = 
0.97 ±0.16). Although rs^i is closer to the average value 
seen in quasar-host galaxies than in SMGs, given that 
similar values are seen in some \ow-z ULIRGs that lack 
AGN, rs^i ~ 1 is not inconsistent with previous observa- 
tional limits on the AGN strength in J14011. J14011 may 
be experiencing a short period of highly- concentrated 
star formation aff ecting the gas excitation (e.g., Yao et al. 



ent is larger than the galaxy wide- aver age, implying that 
J14011's molecular ISM comprises numerous unresolved 
molecular clouds. 

Based on the agreement between the metallicity- 
dependent ver sion of the CO-to-H2 conversion factor 
([2012 ) and the canonical "disk" 
LVG analysis, we use aco = 
^ to derive a molecul ar gas mass 



from Narayanan et al. 
favored" 



value favored by our 
4.6M0(Kkms-ipc2)- 
of (1. 9±0.3) X 10^^ Mq (corr ected for /i 
2005|). Using Ha data from Tecza et al 
W 



3.5: 



(2004) 



S mall et al 
corrected 



2003 lono et al.| |2009), indicating that rs^i ^ 1 is not a 
pertect discriminant between starbursts and quasar host 
galaxies at high redshift. We find that the position of 
the CO (1-0) emission agrees with those of pre viously 
detected mid- J lines (Downes & Solomon 2003). The 
CO (1-0) emission is spatially extended from the south- 
east to northwest, extending between the Jls and Jin 
components previously identified at optical wavelengths, 
and showing two peaks that are slightly offset in velocity. 
A face-on interpretation is favored by our LVG models 
and is consistent with J14011's narrow line width. 

Although the results of the LVG modeling are highly 
degenerate and we are restricted to a single-phase anal- 
ysis, our determination of the gas kinetic temperature 
is consistent with the gas temperature measured using 
C I (Walter et al. 2011) and with the dust tempera- 
ture (Wu et al.„ 2009). Our best-fit models give nu^ 



the SFR measured in the I^IK, we are able to com- 
pare the spatially-resolved Schmidt-Kennicutt relation in 
J14011 to those of other galaxies. We find that J14011 
falls above the star formation rate vs. molecular gas sur- 
face density relation seen in normal star-forming galax- 
ies (e.g., Bigiel et al.||2008|), even when we use the local 
"disk" value of the CO-to-H2 conversion factor. This re- 
sult indicates that the observed offset of starburst galax- 
ies from the local universe Schmidt-Kennicutt law is not 
solely due to the use of tw o different gas mass conversion 



factors, in agreement with Genzel et al. (2010). J14011 is 



17.00±0.25 



cm 



' km 



lO'^-lO^ cm-"^ and Nco/Av = 10 
These densities only produce reasonable [CO/H2] values 
and satisfy dynamical constraints if the velocity gradi- 



consistent in its star formation relation with other SMGs 
and starbursts when we respect different authors' choices 
of gas mass conversion factor. The near-unity index of 
the power law fit to J14011's points is comparable to 
indices measured in the local universe. 
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